Double resonances in Borromean heteronuclear triatomic systems 
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We investigate the occurrence of Borromean three-body continuum s-wave resonances, in an aa/S 
system for large negative two-body scattering lengths. The energy and width are determined by 
a scaUng function with arguments given by energy ratios of the two-body virtual state subsystem 
energies with the shallowest three-body bound state. The Borromean continuum resonances emerg- 
ing from Efimov states present a peculiar behavior for trapped ultracold atoms near a Feshbach 
resonance: two resonances with equal energies at different values of the scattering length. The 
corresponding three-body recombination peaks should merge as the temperature is raised, with one 
moving towards lower values of the scattering length as the other moves to larger values. 

PACS numbers: 21.45.-v, 36.40.-c, 34.10. -fx, 03.75.-b 



The search for Efimov states was favored in recent 
years with the development of new techniques in cold- 
atom laboratories, such as laser-cooling mechanisms and 
the control of two-body interactions by Feshbach reso- 
nance techniques (for recent reports, see [l|, \^ and ref- 
erences therein). At the Innsbruck laboratory, Efimov 
states are identified in scattering using ultracold gas of 
Cesium atoms [3| at 10 nanokelvin. And, more recently, 
at the laboratory of Universita di Firenze [4!] , by using a 
mixture of ^^K and *^Rb it was identified the formation 
of Efimov resonances in the two heteroatomic channels 
for KKRb and KRbRb. The asymmetry of masses may 
turn more favorable the observation of Efimov states in 
comparison to the case of three identical bosons (see e.g. 



Generic systems of the type aa/3 can be classified 
according to the interaction of their two-body subsys- 
tems [a, y. In particular, Borromean systems have been 
studied in |9| . Progress on calcula ting three-body contin- 
uum states and resonance decay jl2| are impressive and 
they give a step further in our precise understanding of 
these states. Following these novel methods, applied so 
far in the nuclear physics context (e.g. the decay of 0+ 
state of ^^C in three a's), a picture of resonances aris- 
ing from Efimov states are highly demanded in view of 
the tremendous advances on the experimental setups of 
ultracold atoms with tunable interactions. Some proper- 
ties of these states under the change in the interaction 
parameters are well known, e.g., when at least one of the 
two-body subsystems is bound, an Efimov state emerges 
from the two-body analytical cut by changing the two- 
body binding energy [1^. Such state, in the second en- 
ergy sheet, is given by a pole in the real axis and identified 
as a virtual state. In the case of a Borromean three-boson 
system there is no two-body cut, and the Efimov state 
arrives from a continuum three-body resonance [11|. As 
systems of ultracold atoms with different species offer 



larger opportunities for universal physics than the equal 
boson case, it is timely to study the fate of Borromean 
states in heteronuclear systems with tunable interactions. 

In the present work, we extend the analysis of Ref. jll|. 
to Borromean aaf3 systems (negative scattering lengths) 
with two kind of particles. It is shown how the Efimov 
states emerge from three-body s-wave continuum reso- 
nances. As the absolute values of the scattering lengths 
(given, respectively, by \aaa\ and/or lac/jj) increase, an 
existing three-body continuum resonance disappears at 
the scattering threshold with the formation of an Efimov 
bound state. This process can be replicated by further 
increasing one or both two-body scattering lengths, with 
the formation of more Efimov excited states. In the ex- 
act Efimov limit, a tower with infinite excited states is 
produced when both, \aaa\ and |a„^|, turn to be infinite. 

The parametric region, for which exists a three-body 
resonance, is presented in the form of a scaling function 
having scattering lengths measured in units of a three- 

r n -1 

body length identified with 



\fm^JhF^ 



where rtic^ 



is the mass of the a particle and i?3 is the three-body 
bound-state that we use as our three-body scaling param- 
eter. In the following, our units will be such that fi = \ 
and rria = 1, with the mass ratio defined as A = m^/niQ. 
In general terms, the quantum description of such large 
and weakly bound systems is universal and can be defined 
by few physical scales, despite the range and details of the 
pairwise interaction. All the detailed information about 
the short-range force, beyond the low-energy two-body 
observables, is retained in only one three-body physical 
information in the limit of zero-range interaction [5, 6| . 
The two-body scales are defined by the energies of the 
two-body virtual states, Eaa and Eap, or by the cor- 
respondingly scattering lengths Oqq ~ —1/ ^J\Eaa\ and 
aa[j ~ —1/ ^J\Eap\. Such two-body scales are given in 

/^, with i?3 de- 



-ll,/\E^\ 
units of a three-body length scale l/^ 



fined as the binding energy of the shallowest three-body 
state. In the present case, the scaling function for the 
continuum resonance energy E^ can be written as 



E-i — B-i £ 



(1) 



The 



corresponding critical boundary, given by 
•\/S3, Oc^-s/Sij; a) = 0, in the parametric plane 

lEctal/B^, ^yiEapl/B^ j , was numerically verified in 

Ref. [l3|, where all possibilities for bound and virtual 
two-body states were considered. We note that the 
calculations done throughout this work were performed 
by considering only two-body virtual states (negative 
scattering lengths). Once the critical boundary is 
crossed in the direction of infinity negative scattering 
lengths, an existing continuum three-body resonance 
will move to an Efimov bound state. Such continuum 
three-body state was indeed observed in the particular 
case of identical bosons, through the resonant behavior 
of three-body recombination rate in an ultracold trapped 
gas of Caesium [Sj] . 

The presence of Efimov resonances has been also re- 
ported in Ref. 4] in an ultracold mixture of "^^K and 
*^Rb gases, through the three-body collision of KKRb 
and RbRbK atoms. In these cases, two resonantly inter- 
acting pairs for positive or negative scattering lengths are 
sufficient to allow Efimov states. Such experiment sup- 
ports the prediction made in Ref. |13| . which is resumed 
in a critical boundary in the parametric plane defined 
by (±^/\Eaa\/ B3, ±^/\Ea/3\/ Bsj , for the existence of at 

least one Efimov bound state. Three-body heteronuclear 
mixtures with resonant interspecies were also been stud- 
ied in |14| using zero-range interaction, where the critical 
conditions for formation of Efimov states at threshold are 
discussed. The case presented in this work of continuum 
resonances were left open. 

In order to calculate the scaling function ([1]), we use 
subtracted Faddeev equations with zero-range interac- 
tions for a three-body system composed by two identical 
bosons, a, and a third one, l3. The following equations, 
presented for continuum resonant states, were already 
writteiL in_a similar form in the context of halo nuclei sys- 
tems [Sl, |lO|. After partial- wave projection, the s— wave 
coupled integral equations of the Faddeev spectator func- 
tions Xaa and Xa/3j can be written as 

Xaa{q)=^TVTaa{q;E3) p dp dzGi{q, p, z; E3)Xaf3{p) , 

/■oc pi 

Xcii{q)=2TTTc,fi{q;Ez) j p^dp dz x (2) 

X [Gi(p,<j, 2;£;3)Xqq(p) +G2{q,P,z;E3)xai3{p)] , 

where the two-body amplitudes, Taa and Tap, and the 
subtracted Green functions, Gi and G2, are defined in 
Ref. jlO|. The subtraction energy is of no physical im- 
portance and it is let to 00. 




1 10 

10' [|E l/Sf"] 



FIG. 1: Real (upper frame) and imaginary (lower frame) parts 
of the three-body complex energy as a function of the two- 
body virtual state energy with E^a = -Eq/s- The results, 
obtained for A*' = 1, are given for mass ratios A —1, 9, 200. 

The resonance energies are calculated from the cou- 
pled eqs. ^ using a contour deformation method in the 
complex momentum plane [15| . The momentum variable 
as p appears as a function of the deformation angle 9 and 
written as p = |p|e~'^, with < 6 < 7r/4. The second 
energy sheet is revealed when the momentum integra- 
tion path is deformed by a contour along the real axis 
to the complex plane with a fixed rotation angle chosen 
to place the contour path far from the scattering sin- 
gularities. That exposes the complex resonance energy 
pole of the scattering matrix. For large enough 9, the 
solution of the coupled eqs. 1^ in the complex energy 
plane is found for tan(26') > —Im{E-i)/Re{E^), where 
Re{Ej^) represents the resonance energy and Im{E3) its 
half width. 

In the limit when the two-body energies are equal to 
zero, an infinite number of Efimov states emerges from 
the solution of the coupled eqs. 1^. At this point, 
by varying the two-body energies, a given energy i?3 
will have zero imaginary part, becoming purely real and 
negative. The N*^ Efimov bound state is defined by 
E^ ' = —B^ ' with A^ = indicating the ground state. 
When the A^*'' Efimov bound state is moved to a reso- 
nance its complex energy will be denoted hy E^ . 

The solution of the homogeneous coupled eqs. (I2|) in 
the complex rotated contour for complex energies gives 
the position and width of the continuum three-body res- 
onances, allowing to construct numerically the scaling 
function ([ij. In practice, we note that one can ap- 
proach quite fast the scaling limit |16j or the limit cy- 
cle [17| for N ^ 00, such that we will present results 
for A^ = 1. In two frames, we show in Fig. [T] the real 



and imaginary parts of the complex energy E. 
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FIG. 2: Given the TV— th three-body resonant energy in units 
of the (A'^ — 1)— th bound state energy B^ ~ , it is plotted 



~Im{Er>)/B., 



,(iV-l) 



function of \J Re{Ef'^)/B. 



for A'^ = 1, Eaa = Eafi and mass-ratios A 



(N-l) 
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1, 9 and 200. 



functions of the two-body binding energy Eaa, for the 
case that Eaa — Eap- AU such energy quantities are 
dimensionless, given in units of the bound-state energy 



B 



(N-l) 



So, in Fig.[Ha), we have Re{E^ 



(N) 



)/B: 



{N-l) 



sus Eaa/Bl,"-'>; and, in Fig. [Hb), /m(sf ))/sf "^^ 
versus Eaa/B^ 



(N-l). 

. In Fig. [ija) , in the positive part of 
the pfot we have the TV— th resonance energy given by 



..u^xv^^^^)/S^^"^\ When the ratio Eaa/B^f'^^ .o ^^- 
creased, this resonance turns into a bound state and the 
negative part of the pfot gives the ratio of the energies of 
two consecutive states (Bg ') and {B^ ~ ') three-body 
bound state for each system. This transition from a res- 
onant to a bound state folfows the same general behav- 
ior as verified in the case of three identical bosons |ll| . 

at 



For ^ = 1, 9 and 200, the values of Eaa/B^/' ^^ 
which we have transitions from resonances to bound- 
states are, respectively, given by 0.876xl0~'^, 1.67xl0~'^ 
and 2.15x10"'^. These values correspond to the ones 
given the position of the boundary curve in the para- 
metric space (—\/\E^^^^^J/B3, —y^lEa^l/B^j (for the case 

Eaa = Eai3) [13| that Separates the excited bound Efi- 
mov state from the continuum resonance. 

In Fig. [TJ our plots are determined only for iV = 1, 
considering that results for higher A'^'s should coincide 
with the ones obtained for A^ = 1. The energy value of 
the resonance (real part) grows, as the virtual two-body 
energies increases (by reducing the absolute values of the 
scattering lengths), up to a maximum and then decreases, 
while the imaginary part always increases. This makes 
the resonance to dive deeper in the second energy sheet, 
as it is clearly shown in Fig. ([2]). 

The results shown in Fig.[l]suggest that in an ultracold 
mixture of heteronuclear atoms the increase of temper- 
ature moves a Borromean Efimov continuum resonance 
towards smaller absolute values of the scattering length 
(larger two-body virtual energies), while the width [Imag- 



inary part of -E3 ] increases. Indeed, an experiment with 
ultracold Caesium atoms [3J indicated that the recombi- 
nation peak due to the continuum triatomic resonance 
moves towards smaller values of \a\ as temperature is 
raised, therefore such resonance is identified with a state 
that just arises from an Efimov bound one that dives 
into the three-boson continuum. Such behavior is a clear 
consequence of the path followed by the real part of the 
resonance energy that increases as \a\ is decreased (see 



Fig.IUandref. [18[) 



The maximum value reached by the real part of the 
three-body energy, the position of the resonance, pre- 
sumably creates a curious effect that one could observe 
in the three-body recombination. As the position of the 
resonance returns to zero {Re{E^ ^) — > 0), it may be pos- 
sible to detect, at a given temperature, the occurrence of 
two peaks in the three-body recombination, correspond- 
ing to two values of the scattering length. The wider 
peak should be at smaller values of the absolute value of 
the scattering length, as one can easily verify by looking 
at Fig. [B 
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FIG. 3: Boundari es for the exist ence of a resonance as a scal- 
ing function of ylEa/il/B,^ ~ versus \/\Eaa\/B.j ~ . 



In Fig. 13] we show a scaling function for E. 
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of y\Eai3\/B^3 versus y\Eaa\/B^^ ', i.e., these 

curves are the boundaries for the existence of a resonance, 
as already obtained in Ref.[l3J. In the lower section of 
the figure, we have three plots carrying the points where 
a resonance begins to appear, Re{E^ ^) = 0, in Fig. [1] 
and, in the upper section of the figure, we have the cor- 
responding curves carrying the points where the origin is 
reached with a finite value of the width. The boundary 
curves indicate that, at ultralow temperatures, a mix- 
ture of heteronuclear gas with two species can have two 
three-body recombination peaks associated with one Efi- 
mov state. As aap is modified with aaa kept fixed, a 
wider resonance appears when the upper curve in Fig. [3] 



is crossed, that happens for sniaher values of \aap\\ with 
a thinner one occurring for larger values. This last case 
corresponds to an excited bound Efimov state diving into 
the continuum. As T is raised, the two peaks tend to 
merge, as the maximum of the resonance is approached. 

We found two continuum resonances for Borromean 
systems at different values of the scattering length and 
same position in energy and different widths, for het- 
eronuclear and homonuclear systems. The absolute max- 
imum for the resonant energy (associated with the sys- 
tem temperature) occurs in both cases. As we observe 
in Fig. [1] the resonance energy attains a maximum value 
and then decreases as \a\ is moved further towards small 
values. Therefore, two resonances are located at the 
same energy for different values of the negative scatter- 
ing length. It is expected that the corresponding pair of 
three-body recombination peaks merges as the temper- 
ature is raised, one peak moves towards lower values of 
the scattering length while the other one to larger val- 
ues. Using the numerical values of the scattering length 
for homonuclear systems, one gets that the ratio between 
the scattering lengths where the real part of the reso- 
nance energy is identical stays between 1 and 0.25 as 
one easily gets by inspection of Fig. [T] Therefore, a sec- 
ond three-body recombination peak appears in a region 
where the four-body recombination peak is present with 
a scattering length ratio of 0.43 [2,ll9|, posing an experi- 
mental challenge for separating the three- and four-body 
recombination peaks. 

In conclusion, we studied the universal properties 
of three-body continuum resonances for heteronuclear 
three-particle system with two species of atoms, in a Bor- 
romean configuration such that all the two-body subsys- 
tems are unbound, i.e., each pair has a virtual state close 
to the two-body scattering threshold. It is shown how 
a Borromean Efimov excited bound state turns out to a 
resonant state by tuning the virtual two-body subsystem 
energies or scattering lengths, with all energies written 
in units of the next deeper shallowest Efimov state en- 
ergy. The resonance position and width for the decay 
into the continuum are obtained as universal scaling func- 
tions (limit cycle) of the dimensionless ratios of the two 
and three-body scales, which are calculated numerically 
within a zero-range renormalized three-body model. We 
discussed how the continuum resonances for heteronu- 
clear Borromean systems may be observed in mixtures 
of ultracold gases. The fingerprint of the observation of 
the scaling considered in the present work, for a Bor- 
romean system, is that such resonances will be given by 
their characteristic displacement with temperature T: as 
T raises the recombination peak, that comes from an 
Efimov state diving into the continuum, should move to 
lower absolute values of the scattering lengths with an 
increasing width. Moreover, we found theoretically two 
recombination peaks at T = 0, with different widths and 
positioned at different values of the scattering lengths. 



For a fixed aa scattering length, with the assumption 
that three-body scale does not move, the wider resonance 
should appear for smaller absolute values of the a/3 scat- 
tering length. As T is raised the two recombination peaks 
tend to merge. This physical effect may be observed in 
ultracold trapped heteronuclear systems near a Feshbach 
resonance, where the properties of the continuum reso- 
nances offer a rich structure to be explored by tuning 
the large and negative scattering lengths. It is an excit- 
ing possibility to be probed in actual experiments with 
mixtures of ultracold atomic gases. 
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